ON 



o 
O 

On 
(N 

O 

o 

Or 



^3 

o 
o 



> 

m 
o 

On 
-i— > 



i 

c 

o 
o 



X 



Magnetic impurity induced id xy component and spontaneous time reversal and parity 
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It is shown that the patches of complex d xy component are generated around magnetic impurity 
in the presence of the coupling between orbital moment of the condensate and impurity spin S z . 
Locally induced d xy gap leads to the fully gapped quasiparticle spectrum near impurity. The phase 
of the induced +id xy component is determined by impurity spin and averages to zero at high 
temperatures. It is suggested that at low temperature the well defined patches of d xy are formed 
and they can undergo a phase transition into phase locked state via Josephson effect. Violation 
of time-reversal symmetry and parity occurs spontaneously via the second order transition. In the 
ordered phase both impurity magnetization and d xy component of the order parameter develop and 
are proportional to each other. 

PACS numbers: 74.25.Bt, 74.20De, 74.62.Dh 



It is well known that magnetic impurities destroy 
the singlet superconducting state due to spin scattering 
which breaks pair singlets M . In the case of the gapless 
(with the nodes of the gap) d-wave superconductor both 
magnetic and nonmagnetic impurities produce a finite 
density of states at zero energy. These effects are a simple 
and direct consequence of the lifetime effects produced by 
impurities. These are well known "incoherent" effects of 
impurities in unconventional superconductors. After re- 
cent experiments by Movshovich et.al. we are led to 
believe that another phenomenon is possible, namely the 
transition to the second superconducting phase as a re- 
sult of condensate interactions with magnetic impurities. 
The time reversal violating state is formed at low energy 
and the order parameter of the new phase is d x 2_ v 2 + id xy 
(d + id). In this phase the impurity spins acquire nonzero 
spin density along z-axis, i.e. out of plane. 

The physical origin of the instability comes from the 
fact that the d + id state has an orbital moment which 
couples to the magnetic impurity spins. The relevant 
interaction is the L Z S Z coupling between impurity spin 
S z and the conduction electron orbital moment L z : 



Hi, 



fJ 



d r 



(1) 



where g - is the coupling constant, ip ra - the electron 
annihilation operator and summation is over impurity 
sites i. In the pure phase one can think of d-wave state 
as an equal admixture of the orbital moment L z = ±2 
pairs: 



A (6) = A cos2O 



Ac 



exp(2i0) + exp(-2i0)) (2) 



Here O is the 2D planar angle of the momentum on the 
Fermi surface, Ao is the magnitude of the d x 2_ y 2 compo- 
nent. We consider 2D d x 2_ y 2 superconductor, motivated 



by the layered structure of the cuprates. In the presence 
of the (ferromagnetically) ordered impurity spins S z the 
coefficients of the L z = ±2 components will shift linearly 
in S z with opposite signs: 



A (6) 



A 



((l+ ff 5 z )exp(2ie) 



+(l-ff£ a )««p(-2ie)) = Ao(e) + iS,Ai(e) (3) 

where Ai(O) cx g/2 sin 20 - is the d xy component. The 
relative phase 7r/2 of these two order parameters comes 
out naturally because d + id state has a noncompensated 
orbital moment L z = +2. 

Here I argue that time reversal (T) and parity (P) 
symmetries can be broken spontaneously in the bulk of 
the d-wave state due to coupling to the impurity spins. 
Original d x 2_ y 2 is unstable towards the formation of the 
bulk d x 2_ y 2 + id xy phase. In the new phase both the 
spontaneous magnetization of impurity spins and sec- 
ond component of the order parameter are developed 
simultaneously. To show how complex d xy component 
appears, I first consider the single magnetic impurity 
and find that the spin-orbit interaction between impu- 
rity spin and orbital moment of the conduction electron 
generate a finite complex d xy anomalous amplitude near 
impurity in d x 2_ y 2 state. This patch of d xy state is 
formed near impurity site, as long as d x 2_ y i amplitude 
is finite, and has a spatial extend of coherence length 
£o = 20A. It is therefore possible for these patches to 
form a long range phase coherent state at some lower 
temperature as a result of Josephson tunneling between 
different patches. I also present a macroscopic Ginzburg- 
Landau functional (GL) and find that there is a linear 
coupling between the original d x 2_ y 2 order parameter 
A (6) = A cos26 and the spontaneously induced d xy 
component: Ai(0) = Ai sin 20. The GL functional con- 
tains the linear coupling term: 
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Fint — — — (A Ai 



h.c.)S z 



(4) 



where b oc ni mp g is the macroscopic coupling constant, 
n-imp is the impurity concentration per unit cell of lin- 
ear size a, Ao,i, g/a have dimension of energy. The time 
reversal violation is natural in this case as it allows the 
order parameter Ao + iAi to couple directly to the impu- 
rity spin. This coupling is possible only for d+id and not 
for d+is symmetry of the order parameter. From the GL 
description it follows that instability develops as a sec- 
ond order phase transition where both the out-of-plane 
magnetization S z and d xy component developed together 
and are proportional to each other ||. 

Recent experimental observation of the surface- 
induced time-reversal violating state in YBCO suggests 
that the secondary component of the order parameter 
(d+is) can be induced Theoretical explanation, based 
on surface- induced Andreev states has been suggested by 
Sauls and co-workers ||. The source of the secondary 
component is the bending of the original d x 2_ y 2 order 
parameter at the surface. 

In a different approach Laughlin || argued that the 
d x 2_ y 2 state is unstable towards d + id state in the bulk 
in the perpendicular magnetic field at low enough tem- 
peratures. The time reversal and parity are broken by 
external field in this case. The linear coupling of the sec- 
ondary order parameter to the external field is central 
to his consideration and results in the first order phase 
transition into d+id state. This transition was suggested 
to be responsible for the kink-like feature in the thermal 
conductivity in experiments by Krishana et.al. |Q 

Recent experiments reported the anomaly in the ther- 
mal conductivity in Bi2212 at low temperatures: the 
thermal conductivity of the Bi2212 with Ni impurities 
was observed to have a sharp reduction at T* = 200mK 
Pj . These experimental data indicate the possible super- 
conducting phase transition in the Bi2212 in the pres- 
ence of the magnetic impurities in a certain concentra- 
tion range. So far the transition has been seen only in 
the samples with magnetic impurities, e.g. Ni as op- 
posed to the nonmagnetic impurities such as Zn It 
was reported that the feature in the thermal conduc- 
tivity is completely suppressed by applying the field of 
H ~ 200Gauss. The low field and the fact that feature 
disappears is consistent with the superconducting transi- 
tion into second phase. Results presented here might be 
relevant for the experimentally observed transition at T* 
in Bi2212 with Ni. 

1. Single magnetic impurity and d xy patch. 

I begin by considering one impurity at site r, = and 
interacting with conduction electrons via Hi nt in Eq. ([!]) . 
Similar to the approach of Q , one can find the anomalous 
propagator in the presence of the single impurity scatter- 
ing potential: F WB (k, k') = F°Jk)5(k- k') + F^Jk, k'), 



are the pure system propagators, (k, k') is the correc- 
tion due to impurity scattering, k = (fc, 6) are the mag- 
nitude and angle of the momentum k on the cylindrical 
Fermi surface, u n is Matsubara frequency and £k = e^ — fi 
is the quasiparticle energy, counted form the Fermi sur- 
face. We take S to be a classical variable and ignore spin 
flip scattering. To linear order in small g one finds : 

i£,(k,k') = -^2 7 ^ 3 5 ^ G°, l (k)F« Il (k')fc^ (5) 

Where F* (k, k') is the function of incoming and outgo- 
ing momenta because of broken translational symmetry. 
Upon integrating F 1 over k' and going to integrated over 
£k propagator one finds: 

<,(©)= / N dfi k FZ, n (k)=iA uln (N gS z )(N A ) S m2Q (6) 



Here A Wn ~ k F ^hx{- 



w 



cos- 



26) t 



is the model dependent coupling, W is the energy cutoff, 
No is the Density of States at the Fermi surface and ()e 
stands for Fermi surface averaging. 

The angle dependence of F* (8) ~ igS z sin 20 cx k x k y 
is the one of d xy . Together with the fact that this am- 
plitude is complex it indicates the existence of the id xy 
component in the vicinity of magnetic impurity. This 
result also shows that incoming d x 2_ y 2-w&ve state elec- 
trons have a finite amplitude, linear in gS z , to be scat- 
tered into the d xy outgoing state via the L Z S Z coupling. 
From the solution Eq. (|]) it is easy to see that the typical 
size of the patch, ignoring nodal directions, is given by 
superconducting coherence length £o = 20A For the rel- 
evant concentration of Ni ni mp ~ 1% the Ni-Ni distance 
is about 35A. The patches are thus well overlapping in 
this limit making phase ordering due to tunneling form 
patch to patch possible. These patches work as a micro- 
scopic seed of d xy component which grows into true long 
range state at low temperatures T < T*. Similar result 
for the d xy patches in the mixed state of pure Bi2212 was 
shown in ||, where the role of impurities is assumed by 
vortices. 

Important observable consequence of the local d xy com- 
ponent near impurity is that the gap, as seen in STM 
tunneling near Ni impurity, will increase and the low en- 
ergy part of Density of States will be suppressed because 
of finite gap everywhere on the Fermi surface, as opposed 
to nodes for pure d x 2_ y 2 state. The increase of the gap 
in the Ni-doped Bi2212, compared to the pure case, was 
observed in STM tunneling @ 

Next, I consider simple example of magnetically or- 
dered state: the ferromagnetically ordered impurity 
spins. This does not have to be the case experimentally 
but makes the point clearer. Calculation, similar to the 
one above, yields: 



where F° 



An cos 20 



"J+5k+A2cos 2 2e 



G° 



w2-K k +A= cos 2 2G 
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gS z G° (k)[9kXk] z F° (k) (7) 
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The existence of the homogeneous d xy component 
F,J (k) ~ — iriimp^-gSzAo sin 28 is evident from this 
equation. The relative phase of d xy component with re- 
spect to d x 2_ y 2 is determined by the sign of S z . The gap 
Ai has to be determined seflconsistently provided there is 
interaction in xy channel. This interaction does not have 
to be attractive, since d x 2_ y 2 plays the role of the source 
and Ai will be generated for any sign of interaction. I 
will assume there is such interaction and results below 
will be expressed in terms of the induced gap Ai. With 
the help of this equation I find for the free energy change 
due to A x : 5F = l/2T£ Wn k i^(k)AJ(k) + h.c. 

SF = -i/2(NoA{)(N Ao)^(gS z )n imp + h.c. (8) 

Eq. (||) together with single impurity result Eq. (^J) are the 
main results of this section. From this equation I find the 
linear term Eq.(^|) with 

b = N 09 n irnp (9) 




FIG. 1. a) Impurity sites with random spins at high tem- 
peratures are shown. Spin flips lead to d xy component aver- 
aged to zero, b) Upon slowing down and freezing of impurity 
spins the patches of d ± id with L z = ±2 state near each 
Ni site are formed, shown with right (left) circulating current 
near each impurity site. At low temperatures Josephson tun- 
neling locks the phase between patches, leading to the global 
d + id state. 

So far the spin flips were ignored. The relative phase 
of second component is determined by the sign of S z . 
At high temperatures T> T c * , when spins are strongly 
fluctuating, the relative phase of the d xy component is 
fluctuating strongly as well. This phenomenon is an in- 
teresting new realization of the superconducting phase 
(d xy ) coupled to the heat bath (fluctuating spins). 

If and when the impurity spins are slowing down or 
even are freezing out then the phase scattering time be- 
comes large and the phase ordering of the patches is pos- 
sible, see Fig.l. Measurements indicate that spin flips of 



Ni spins are slowing down at low temperatures T < 2K . 
Specific heat measurements on Ni- doped Bi2212 indicate 
additional entropy, compared to undoped Bi2212, on the 
order of rii mp R log 3, accumulated near IK |y]]. This 
additional specific heat has a broad maximum around 
IK and linear slope at lower temperatures. The general 
shape of the specific heat, associated with the impurity 
spins is strikingly similar to the specific heat, observed in 
spin glasses pi. Broad peak in the specific heat might 
indicate the glassy behavior of spins at lower tempera- 
tures. 

2. Mean field formulation. 

Below I will ignore the fluctuations in the magnetic 
subsystem and consider simple mean field theory of the 
coupled magnetic impurities and superconducting con- 
densate. From the specific heat measurements we know 
that some type of order (spin-glass or other) might oc- 
cur at T m ~ IK. I will assume that impurities develop 
a ferromagnetic order at some temperature T m . This is 
a drastic oversimplification because of the possible spin- 
glass ordering discussed above. Nevertheless the model 
presented below is useful to understand how the coupling 
between impurity spins and condensate leads to the d+id 
instability of the original state. 

I will consider the GL theory of the secondary super- 
conducting transition: Ao — > Ao + Ai at T*, where 
both Ao,! are homogeneous variables corresponding to 
the macroscopic ordering. The relative phase of Ai with 
respect to the phase of Ao is not fixed and will be de- 
termined by the free energy minimization. Assume that 
the second transition, if at all, occurs at T* <C T c , where 
T c ~ 90K is the first transition temperature. Hence the 
order parameter Ao, which can be assumed to be real, is 
robust and its free energy F(Aq) can not be expanded in 
A . 

The relevant fields to enter the GL functional are: Ao, 
Ai and S z . 

Assuming expansion in powers of small S z ,Ai near 
second transition, the GL functional F — F sc + F magn + 
F int is: 

F sc = F(A ) + c^AxI 2 + a 2 /4|A!| 4 + 

^IViAxlVx^O 
F magn = ^\S Z \ 2 + f\S z \ 4 + f |VA| 2 

F int = -±(A£A 1 -h.c.)S z (10) 

Ao should enter in Fi nt for it to be invariant under the 
global U(l) symmetry A ,i — > A 0j i exp(i9) Homo- 
geneous solution will have lowest energy and therefore 
gradient terms are take to be zero hereafter. 

All but Fi n t terms in the free energy Eq.(|lO|) are pos- 
itive and can not produce the instability of the original 
d x 2_ y 2 state. Fi n t can be negative since it is linear in 
Ai and S z and this term is crucial in producing second 
transition. 
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Magnetic energy F magn has a temperature dependent 
coefBcient 

oi(T) = ain imp (T - T m ) (11) 

and vanishes at T m , a\ is dimensionless. With this choice 
I will assume that Ni impurities would order ferromag- 
netically at T m ~ IX in the absence of the interaction 
with condensate. 

Consider F sc . The second and third terms in F 8C de- 
scribe be the energy cost of opening the fully gapped state 
with Ai when interaction prefers to keep node, i.e. pure 
d x 2_ y 2 state. The change in free energy due to secondary 
order parameter is given by the difference in energy of 
quasiparticles before and after Ai component is gener- 
ated. One can calculate the change in the energy of the 
superconductor subjected to the homogeneous external 
d xy source field: H xy = Ax(k)^ k(T ^* k _ CT , where 

< k < 1 is the integration constant. Using standard 
result d K F sc = l/n{H xy ) one finds increase of energy at 
Ai<<T<<A : 

6F SC = ^|A x | 2 + f |Ax| 4 

(12) 

Here a.\ ~ N , and Ai is taken to be constant on the 
Fermi surface, see Eq.(||), p|,|l4|]. 

Fix Ao to be real positive and the relative phase of 
Ai = | Ai| exp(w) without loss of generality. Minimizing 
the functional Eq.(|lO|) I find that the tt/2 relative phase 
of Ai comes out naturally: phase will be determined by 
minimization of energy: 

v = tt/2 sgn{bS z ) (13) 

This choice takes the maximum advantage of the L Z S Z 
coupling and minimization requires a complex order pa- 
rameter Ao + iAi in the low temperature phase. T and 
P are violated spontaneously even with T m = 0. The 
minimization yields: 

S * = ni n imp (r-r m ) sin ^A 1 Ax | 

SF^-^A^- \T-T*\ 2 
T* ~T m + ^ 2 {N Q A ) 2 n lmp 9 ^ (14) 

where I used the Eq.(||) in the last line. This is the main 
result of this paper. Solution Eq.(|lJ) indicates that the 
transition is of the second order with the jump in the 
specific heat. I assumed that |Ao(T)| 2 /a\(T — T m ) has 
a linear temperature slope near T*. It follows from the 
solution Eq([l4j) that: 

1) Even if T m = the ordering will occur at T* = 

47r 2 rii mp (iVoAo) 2 ^^. However the softness of the spin 
system near T m enhances the effect and makes T* > T m 



within this mean field approach. Taking the typical val- 
ues for Bi2212 of A = 450A, A F = 1/N = 3000 A, as- 
suming a\ r*> 1 and taking the characteristic value of spin- 
rbit coupling of Ni g/a ~ 320A', I find T* ~ 0.3A. In the 
real system, if the spin-glass freezing occurs, the freezing 
will occur at first for the spins that are well separated. 
This will preclude the Josephson tunneling between the 
patches, as discussed above. Only at lower temperatures, 
when majority of spins are frozen, the tunneling would be 
able to lock in the superconducting phase. Hence the real 
phase ordering will occur at temperatures, substantially 
lower then the mean field estimated T* < T m . 

2) As the function of impurity concentration two effects 
occur simultaneously. First, condensate density |A | 2 de- 
creases. Second, the suppression of Ai due to increased 
impurity scattering will also lower T* . These effects will 
lead eventually to the disappearance of the transition. 
Quick suppression of the transition temperature T* with 
impurity concentration should be expected. 

3) Strong magnetic field parallel to the layers, H 3> 
Hci,ab ~ 1 Gauss in plane, will suppress the second 
phase. In the field Ni spins will be aligned in the lay- 
ers, linear coupling term on Hi„ t will be zero and d xy 
component will vanish. This effect might explain the 
suppression of the second transition by magnetic field 
H c i ab "C H < H c i c ~ 300 Gauss, seen in experiment 

Weak localization of the quasiparticles [15|, in princi- 
ple, can cause the rapid decrease in the thermal conduc- 
tivity. Experimental facts argue against localization in 
the layers for the following reasons: the field parallel to 
the layers suppresses the observed feature, it disappears 
at higher Ni concentration and the specific heat is in- 
creased near transition temperature. Specific heat and 
thermal conductivity in the field parallel to the layers 
will help to determine how relevant localization of quasi- 
particles is to the experimentally observed transition. 

In the superconducting state with nonzero orbital cur- 
rent L z the dominant fraction of the orbital moment is 
"stored" at the edge of the sample, similar to 3 He — A. 
The edge currents in d + id state and their topological 
characteristics was addressed recently in |||. 

To test the proposed state following experiments can 
be done. The driving mechanism for the second phase 
clearly distinguishes between magnetic and nonmagnetic 
impurities, hence more experiments on Bi2212 with non- 
magnetic impurities will be helpful [0. Theory predicts 
the ordering of Ni moments below T*, and one should be 
able to detect magnetization in fiSR experiments or in 
ac susceptibility. The increased superfluid density due to 
second component translates into the change in the pen- 
etration depth below T* which can be detected. These 
and other experiments will help to resolve if the proposed 
mechanism is correct. 

In conclusion, I presented the mechanism for a sec- 
ond order phase transition of original d-wave state into 
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d + id state with spontaneously broken T and P. In the 
ordered phase both impurity spins S z and d xy component 
of the order parameter develop and are proportional to 
each other. The low temperature phase develops mag- 
netic moment both due to magnetic impurities spin and 
because of the finite angular momentum of d + id state. 

I am grateful to R.B. Laughlin, D.H. Lee, A. Leggett, 
R. Movshovich, M. Salkola and G. Volovik for the useful 
discussions. This work was supported by US DOE. 
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